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ABSTRACT 


Under  Contract  AF  33(6571-  14648.  both  routine  operation  of  the 
Cumberland  Plateau  Seiemologieal  Obaarvatory  and  conelderable  applied  r«- 
eearcb  were  conducted.  During  the  laet  project  year,  on-line  real-time  detec¬ 
tion  and  identification  proceeding  ueing  the  CPO  Aux'liery  Proceeeor  wae  im¬ 
plemented  and  evaluated  at  CPO  for  the  purpoee  of  etudying  automatic  detection 
proceeding.  The  CPO  Auxiliary  Proceeeor  compute*  two  claaeee  of  detection 
output*,  the  Fieher  analyei* -of- v  iriance  etatietic  and  the  Wiener  power  etatii- 
tic,  and  one  clae*  of  identification  output,  the  United  Kingdom  technique. 

Theee  detection  output*  were  compared  on-line  againet  a  fixed  eignal  threshold 
level,  providing  a  ecntinuoue  real-time  "ye*-no"  output  for  eignal.  Howeve*. 
the  iixed-thrcehold  detection  level*  were  initially  difficult  to  determine  accur¬ 
ately  and,  once  determined,  it  wa*  found  that  they  were  highly  non-time  *ta- 
t ionary .  Adaptive  threshold  detector*  mcorporeted  into  the  Auxiliery  Pro¬ 
cessor  could  overcome  the  non-time  atationarity  of  the  thre*hr  d  detector*. 
Off-line  applied  research  wa*  performed  to  support  the  on-Une  -eaearch. 
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SECTION  I 


INTRODUCTION 


In  addition  to  routine  operation  of  the  Cumberland  Plateau 
Seismologies]  Observatory,  considerable  applied  research  has  been  con¬ 
ducted  under  Contract  AF33{657)-14fc48  to  advance  the  unde  retarding  of 
array  processing  technology  applicable  to  smalt  seismic  arrays  used  in 
the  nuclear  detection  and  classification  problem.  During  the  last  project 
year,  under  VT/67Q4  on-line  real-time  detection  and  identification  pro¬ 
cessing  was  implemented  and  evaluated  at  CPO  for  the  primary  purpose 
of  studying  automatic  detection  processing. 

On-line  processing  was  implemented  using  the  CPO  Auxiliary 
Processor  {Figure  I-1J  which  computes  two  classes  of  detection  outputs,  the 
Fisher  analysis -of-varlance  statistic  and  the  Wiener  power  statistic,  and 
one  class  of  identification  output,  the  United  Kingdom  technique.  The  de¬ 
tection  outputs  were  compared  on-line  against  a  fixed  signal  threshold  level 
for  automatic  detection.  From  thie  comparison  a  continuous  real-time 
"yes -no"  output  was  provided  for  signal.  This  output,  along  with  the  de¬ 
tection  and  identification  data,  was  recorded  on  Develocorder  film. 

Contained  in  this  report  is  a  brief  description  of  the  digital 
processing  hardware  and  results  of  the  hardware  evaluation.  Conclusions, 
as  well  as  supporting  data  regarding  the  evaluation  of  on-line  detection  and 
identification  processing,  are  presented  Also,  results  from  off -line  sup¬ 
porting  applied  research  are  covered. 


1-1 


•otenoe  ssnrtsss  division 


0 


Qm  lift*  r#al*tiroe  Identification  *rvd  tulpmtUr  detection 
p,  easing  began  at  Cf*0  *n  Dv>  «mb*  r  tffci  iH*  Installation  and  op^r- 

1U0A  oj  t he  A  *Uia  r r  l*»  *nOr  This  unit  Inlirfic*)  willi  the  CPO  MCF 
pr  -  estor  fr.ro  *Htb  ll  d*fi  **  *11  Input  »t|n*l  diU  Tb»  Auitiliary 
I’fMnior  compute*  two  Unit*  d  Kingdom  output*.  >ne  Tuber  analysis  -of- 
%  *  r  U  r.c  *  *!*<■; i i  output  and  four  Wiener  power  detection  output*  A  "  y*  » - 

ao  4>»(ii ion  (i  r  slgral  It  presented  In  real-time  bucd  upon  the  threshold 
level  del*  lion  output  c<-roparieon  AH  data  I*  displayed  us  Develocorder 
dim 

T h*  primary  r*a*on  lor  implementing  the  Auxiliary  Protestor 
wee  to  rtilertf  the  feasibility  end  effectiveness  of  automatic  defsetion  pro- 
trr*|n|  ti  applied  to  the  rfetectioa  problem  for  email  diameter  short  - 

period  eelemte  array*.  Compart  eon  of  Fisher  autletlc  and  Wiener  power 
del*  lion  protesting  over  large  %  autitiee  of  data,  ae  welt  ae  i6m*  claeel- 
flc atlon  work  «  ting  the  UK  proceaeed  data,  wae  planned. 

After  the  Auxiliary  Processor  design,  conn  ruction  and  in¬ 
stallation  phase  from  alay  Ifh4  through  December  f$Mt,  a  performance  and 
ha  rdware  evaluation  phase  wae  conducted  to  answer  questions  regarding 
feasibility*  reliability,  and  effect  I  vane**,  Concurrent  with  iki#  writ,  off¬ 
line  applied  research  directed  toward  esfending  the  present  processing 
tec  hoe  log  y  and  understanding  wae  prrfdrraed.tn  Dallae, 

A.  HARDWARE  EVALUATION 

Results  m*  the  hardware  evaluation  Indicate  that  the  eyetem  le 
highly  reliable.  Mot  one  component  failure  ei  erred  during  the  JO  December 
*fbb  to  10  AprtJ  Ifbt  operating  period.  The  initial  set  p  and  programing  of 
tfe*  syture  le  difficult  ts4  requires  the  empirical  determination  of  many 
parameters  whi  a  must  be  manaaHy  programed,  fhfe  difficulty  is  particul¬ 
arly  true  of  the  Fieber  process  ismi  knowledge  of  the  Fisher  Intermediate 
term  values  IQ,  K1  and  KJ  roust  be  known  in  order  to  determine  the  optimum 
data  truncate  settings.  Several  recommendations  for  modification*  of  possible 
future  systems  are  presented  In  Section  HI 


B.  ON-LINE  PROCESSING  EVALUATION 


1.  Detection  Processing 

Evaluation  of  the  Auxiliary  Processor  as  a  detection  device 
Indicated  considerable  difficulty  exists  with  the  fixed-threshold  «.  Jtection 
levels  for  the  Pisher  and  Wiener  outputs.  Initially  it  *vai»  difficult  to  deter¬ 
mine  accurately  the  desired  threshold  level.  Once  the  levels  were  deter¬ 
mined.  it  was  found  that  they  were  highly  non-  ime  stationary  which  meant 
that  the  false -alarm  rate  could  not  be  fixed.  Attempts  to  adopt  standard 
procedures  to  update  the  threshold  levels  on  a  daily  (and  sometimes  more 
often)  basis  proved  to  be  inadequate.  The  vaiiations  in  threshold  level 
were  significant  enough  that  the  automatic  detection  outputs  were  for  prac¬ 
tical  purposes  us. '.ess. 

Adaptive  threshold  detectors  should  be  incorporated  in  the 
Auxiliary  Processor.  An  adaption  algorithm  is  presented  in  the  appendix 
which  is  suitable  for  incorporation  into  the  existing  system.  It  is  felt  that 
the  Auxiliary  Processor  can  be  a  significant  tool  in  the  automatic  P-wave 
detection  problem  with  this  addition 

2.  Identification  Processing 

The  two  UK  outputs  were  properly  implemented  in  the  Auxiliary 
Processor  but  are  of  limited  use  at  CPO  for  classification.  First,  this  type 
of  computation  is  designed  for  use  on  a  large  diameter  crossarray  (20  km  has 
been  used).  The  CPO  array,  which  is  3.  6  am  in  diameter,  lacks  sufficient 
directional  resolution  and  violates  the  assumption  that  noise  is  uncorrelated 
across  the  array.  Second,  the  two  outputs  ma^  be  programed  for  only  two 
directions,  severely  limiting  the  class  of  signals  which  may  be  studied. 

Third,  classification  work  requires  the  preservation  of  signal  waveform  for 
all  events.  This  causes  a  bast'’  dynamic  range  conflict  with  the  MCF  which 
is  a  detection  device  requiring  adequate  noise  for  coherent  noise  uuppression. 
The  MCF  Auxiliary  Processor  system  is  limited  to  a  12-bit  (66-db)  dynamic 
range  on  input.  Since  on-line  emphasis  was  placed  on  MCF  signal  extraction 
and  detection  processing,  the  class  of  signals  available  for  study  was  highly 
restricted.  Large  si^iais  of  interest  were  clipped  on  input  or  during  inter¬ 
mediate  computations  and  small  signals  of  interest  (e.  g. ,  NTS  shots)  were 
not  detected  on  the  Deveiocorder  display.  An  adequate  library  of  events  for 
study  was  not  collected  because  of  these  limitations. 

Work  in  this  area  was  subsequently  shifted  to  the  MCF  and 
Auxiliary  Processor  evaluation.  The  identification  processing  technique 
warrants  study  if  the  processor  is  installed  at  a  more  suitable  array  loca¬ 
tion  and  sufficient  events  can  be  accumulated. 
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C.  OFF-LINE  APPLIED  RESEARCH 


I 


Off-line  Dallas -based  supporting  research  was  directed 
primarily  toward  determining  parameter  specifications  fo*  the  Auxiliary 
Processor  program  and  toward  investigating  properties  of  the  Fisher  out¬ 
put.  Two  critical  parameters,  the  integration  gate  length  for  the  detection 
outputs  and  the  low-cut  frequency  filter  specification  necessary  for  pre¬ 
filtering  the  Fisher  input  data,  were  determined  from  studies  of  CPO  signal 
properHes.  Effect  of  correlated  noise  on  the  Fisher  output  was  studied 
empirically  in  relation  to  the  low-cut  lilter  specification.  Compared  to  the 
Wiener  outputs,  signal  attenuation  for  the  Fisher  computation  as  a  function 
of  wavenumber  was  determined  to  be  significantly  greater. 

Thf  variable  threshold  problem  for  a  fixed  false-al&rm  rate 
was  investigated  and  found  to  be  related  to  the  RMS  input  noise  level  for 
the  Fisher  output.  This  i  -rease  could  not  be  related  to  the  effect  of  a 
particular  noise  contributor,  but  it  is  rea-onaole  to  assume  that  there  may 
be  a  strong  correlation  to  the  mantle  P-wave  noise  level. 

D.  CONCLUSIONS 

Much  useful  information  for  advancing  the  state-of-the-art 
real-time  automatic  detection  processing  was  gained  by  the  implementation 
and  evaluation  of  the  CPO  Auxiliary  Processor.  Both  the  Wiener  power  and 
Fisher  statistic  computations  appear  to  be  useful  for  detection  purposes, 
but  completely  automatic  detection  should  be  accomplished  with  an  adaptive 
threshold  device.  Such  a  device  can  be  easily  incorporated  into  the  existing 
hardware  with  a  relatively  minor  modification. 

Little  information  was  gained  from  computation  of  the  UK 
technique  on-line  at  CPO.  Thu  adverse  array  properties  and  dynamic  range 
problem  coupled  with  the  sparseness  of  desired  events  severely  limited  the 
necessary  library  of  events  required  for  this  study. 

Off-line  Dallas  supported  research  complimented  the  on-line 
evaluation  and  provided  a  more  comprehensive  understanding  of  on-line  Fisher 
and  Wiener  automatic  detection  processing. 


SECTION  III 


PROCESSOR  IMPLEMENTATION  AND  EVALUATION 


The  purpose  of  this  section  is  to  provide  a  concise  description 
of  the  GPO  Auxiliary  Processor,  present  results  of  the  hardware  evaluation 
conducted  during  the  on-line  operating  period  and  recommend  system  mod¬ 
ifications  based  upon  experience  gained  through  system  operation. 

A.  SYSTEM  DESCRIPTION 


A  detailed  system  description  may  be  found  in  the  Instruction 
Manual,  Operation  and  Maintenance  Auxiliary  Processor  and  the  Specification 
for  Auxiliary  Processor  Modification,  Advanced  Multi-Channel  Filter. 

1.  Description  of  the  Key  System  Operations  • 


a.  Fisher  Process 


The  Fisher  process  is  a  statistical  signal  detection  process 
which  operates  on  the  filtered  outputs  of  each  MCF  input  channel.  Filtered 
outputs  are  formed  by  the  MCF  in  the  beam-steer  process  using  the  following 
equation. 


whe  re 


(1) 


n 


is  the  value  of  the  most  recent  sample  of  the  data  for 
input  channel  a  (where  k  ranges  from  0  to  k) 


through  a^  are  constants 
1+1  is  the  number  of  filter  points 


By  using  the  CPO  MCF,  the  auxiliary  processor  computes  a 
Fisher  output  according  to  the  following  equation. 


where 


* 


■ym 


Ni  end  N2  are  normalisation  con*tant« 

K  l  is  the  number  of  input  channels 

P  is  the  Fisher  history  length 


Kl,  K2,  and  K3  are  the  intermediate  Fisher  terms  which  are  formed  using 
the  following  equations. 
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The  Fisher  output  computed  by  the  processor  is  a  transform 
of  the  true  Fisher  process  which  is  given  by 


Under  highly  coherent  signal  conditions,  the  denominator  of  F  is  0,  an  un¬ 
defined  state  for  the  hardware  as  implerra  nted.  To  avoid  this  condition 
P  is  transformed  as  follows: 


'  jF® 


- 


■  <  *  y. 


>  * 


\;t ,  '  ; 


-  S* 

y 


i 


C 


0 

B 


#r 


f 


\ 


■  ^ 


Fisher  Output 


F  +  B 
F  +  D 


F  +  D/C 


where 


The  effect  of  this  transform  as  a  function  of  N2  is  illustrated 
in  Figure  Ul-1,  N1  sets  the  maximum  value  (routinely  77?g  or  512^,  the 
largest  possible  output  number  as  a  result  of  the  9-bit  output  register). 
Selection  of  N2  is  based  on  the  F -value,  computed  for  the  site  ambient -noise 
field  (at  CPO,  F=*  2.  2).  This  may  be  determined  operationally  by  finding 
one  or  two  processor-average,  ambient-noise  output  values  over  a  pre¬ 
determined  gate  by  using  the  Fisher  threshold  detectors  and  by  reading  the 
F  value  for  the  known  N2. 


b.  United  Kingdom  Process 


The  United  Kingdom  process  performs  the  0-lag  crosscor¬ 
relation  of  two  beam-steer  outputs.  The  auxiliary  processor  forms  two 
such  outputs  using  the  following  equation. 


UK  Output 


'n-(p-l)  n-(p-l) 


i  /  » 
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AUXUIARY  PROCESSOR  OUTPUT 


where 


beam-steer  output  of  the  MCF  processor 

n- (p-l) 

*  k**1  beam-steer  output  of  the  MCF  processor 

n- (p-  1) 

i  4  k 

P  ■  number  of  UK  history  points 


n  *  mos*  recent  beam-steer  output 


n- 1  *  next  most  recent  beam-steer  output,  etc. 


A  second  UK  output  is  generated  using  outputs  1  and  m,  where 
1  4  m  and  j,  k,  I,  and  m  are  not  equal. 

New  1<K  outputs  are  computed  every  50  msec. 


c.  Wiener  Power  Process 


The  auxiliary  processor  forms  four  Wiener  power  outputs 
using  the  following  equation. 


where 


R  S 


’  E  E  (Am) 


r*l  s*l 


n(s-l)  -  (r-l)S 

the  MCF  output  at  time  m,  truncated  to  11  bits  plus  sign 


the  most  recent  MCF  output  (A_^j  is  the  next  most 
recent  MCI  output,  etc. ) 


sample  interval  sample  interval 
frame  time  50  msec 


window  width  in  sec 


sample  interval 
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The  multichannel  power  processor  output  Is  an  R  x  S  x  50-msec 
window  of  an  S  x  50-meae  sample  interval  of  an  MCF  output, 

d.  T  ireehold  Detector* 

Digital  thra ahold  datactor*  ara  providad  in  the  Auxiliary 
Procassor  for  the  Fisher  output  and  tha  four  power  output*.  Tha  threshold 
detector*  provide  a  eeparate  eignal  output  when  a  monitored  output  equal* 
or  exceed*  a  awUch-aelected  level. 

2.  General  Characteri*tic* 

The  deaign  philosophy  and  general  appearance  of  the  Auxiliary 
Processor  were  modeled  after  the  ba.ic  MCF  processor,  The  unit  consist* 
of  an  80-ln.  ,  aingle-bay  rack  containing  three  drawer*:  the  arithmetic 
drawer,  the  output  drawer,  and  the  controller  drawer  (Figure  1-1).  A*  wUh 
the  CPO  MCF,  all  input  and  output  cabling  paaie*  through  the  top  of  the 

cabinet 

The  apare  rack  space  in  the  Auxiliary  Processor  may  be  used 
for  support  equipment  such  as  the  paper  tape  reader  (PTR),  data  control 
modules,  standard  station  timing  unit,  and  the  rack-mounted  oscilloscope. 


a.  Input  Signals 

All  input  data  signals  are  derived  from  the  basic  MCF  pro¬ 
cessor  as  described  in  the  following  listing, 

•  Wiener  Power  Process 

lhis  process  derives  its  input  from  the  MCFt-4  outputs. 
Of  the  24  available  bite  of  daU.  9  are  aelactad  using  of 
the  input-data  truncation  awitch. 

•  Fiaher  Pieces? 

Single -channel  data  are  accepted  from  tha  MCFO  output. 
Thia  allows  single -channel  prefiltaring  of  data  prior 
to  computation  of  tha  Fisher  statistic.  Of  tha  available 
24  bits  of  data,  9  are  selected  through  the  Fisftar  Input- 
data  truncation  switch. 


•  UK  Process 


The  input  is  derived  from  two  selectable  beam-steer 
outputs.  Nine  of  the  total  24  available  bits  are  se¬ 
lected  by  the  UK  input-data  truncation  switch. 

b.  Output  Signals 

The  digital- to-analog  converters  used  in  the  Auxiliary  Pro¬ 
cessor  are  of  the  same  type  as  those  used  in  the  MCF  processor.  For  the 
Fisher,  Wiener  power,  and  UK  outputs,  9  bits  plus  sign  are  selected  from 
the  2S-bit  output  register  using  individual  output  truncation  switches,  thus 
giving  an  analog-type  gain  control  in  6-db  increments. 


The  processor  outputs  are  summarized  as 


e  One  Fisher  output  trace 

*  Two  UK  output  traces  corresponding  to 
two  programed  area  locations 

e  Four  W iener  power  traces  corresponding 
to  the  MCF1-4  output  traces 

W 

*  One  Fisher  threshold  trace  corresponding 
to  the  Fisher  output  trace 

*  Four  Wiener  threshold  traces  corresponding 
to  the  four  Wiener  tower  traces 

e.  Program  Selection 

The  following  key  variable  programing  modes  are  offered 
by  means  of  panel -mounted  switches. 


*  Fisher  Process 

Normalization  constants  Ni  and  N2  are 
variable  from  0  to  777g,  History  length 
Cgate  length  of  computation)  may  be  selected 
from  0  to  999  points. 
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•  UK  ^Process 

Hiatory  length  ia  variable  from  0  to  999 
point* i  and  selection  of  the  4  beam-ateera 
(2  each  for  UKO  and  UK1)  ia  provided. 

•  Wiener  Power  Procaaa 

Hiatory  length  apecifiable  In  R  intervala, 

0  to  99,  and  S  samples,  0  to  99,  where 
the  gate  ia  determined  by  R  intervala  of 
S  aampiea  each. 

•  Threahold  Detectora 

Independently  variable  threahold  levela 
are  programme  from  0  to  777g. 

B.  HARDWARE  EVALUATION 

1.  Reliability 

The  CPO  Auxiliary  Proceaaor  proved  to  be  a  highly  reliable 
digital  proceaaing  device.  During  the  time  the  unit  waa  operated  at  CPO 
from  30  December  1966  to  10  April  1967,  not  one  failure  occurred.  Thia 
represents  a  total  operating  time  of  approximately  2900  ties.  Normally 
the  "infant  mortality"  period,  during  which  most  failures  are  expected 
to  occur  because  of  marginal  components,  etc.  t  la  set  at  around  2000  hra. 

In  future  operations  this  unit  should  continue  to  be  quits 
reliable  since  it  consists  entirely  of  solid-state  hardware  with  the  exception 
of  the  blower  motor. 

2.  Maintainability 

Little  maintenance  experience  was  gained  on  the  unit  du*  to 
the  high  reliability  experience.  Once  installation  and  operational  checkout 
waa  complete,  processor  maintenance  waa  limited  to  routine  operation# 
which  include  periodic  cleaning  of  the  blower  motor  air  filter  and  the  dally 
step  teat  (thia  teat  is  also  a  part  of  the  MCF  maintenance  procedure  and 
completely  verifies  the  Auxiliary  Proceaaor  logic  and  arithimeticj. 


During  the  en-sile  installation  and  operational  checkout  phase, 
iomi  maintenance  knowledge  was  gained.  The  following  empiric*?  con¬ 
cision*  can  be  reached  based  upon  Ihl*  limited  experience: 

•  The  processor  can  o*  maintained  by  an  indi  vi  du* 
trained  in  MCF  maintenance  if  be  read*  and  thoroughly 
understand#  the  Auxiliary  Proc«**or  Handbook. 

•  The  computer  wire  Hat  provided  with  the  *y*t*m  te  *de 
quale,  but  ■omewhat  tumbtriom*  to  uae.  Document* 
covering  the  MCF  modification  and  interface  were  pro¬ 
vided  a*  "red-lined"  change*  on  the  original  MCF  doc 
umentetion  and,  where  applicable,  reference*  to  this 
change  ll*t  indicated, 

J.  Operation 

Once  the  MCF  Auxiliary  Proce»*or  eyetem  wa»  initially 
programed,  adjustment  wae  limited  to  daily  change*  of  the  threshold-level 
df tiction  twttchti  (or  the1  Fiihf  f  *nd  Witn^r  output*# 

AH  MCF^AuxiHary  Proctiio r  oj/trations  w«rt  ptrformad  by 
station  analyst*  and  technician*  who  had  been  briefed  on  the  *y*tem  operator 
control*. 

4  Progr  ming 

Considerable  difficulty  wa*  experienced  in  initially  programing 
tb*  Auxiliary  Processor  system.  Some  eauae*  for  this  difficulty  are  a* 
follows? 

#  Determination  of  the  input  data  truncation  *witeh 
setting  for  the  Fisher  Procea*  requires  kot  sitrfj* 
at  KI.  K2.  and  KJ,  Thia  information  i«  not  readily 
available, 

•  Determination  of  the  Fisher  intermediate  summation 
truncate  switch  setting  al*o  requires  knowledge  of  the 

•  is*  of  the  intermediate  terms  KI,  Kf  and  KJ> 

»  Determination  of  the  Fish* r  traa*forra  constant  H2 
initially  required  detailed  knowledge  Of  the  F  values 
far  the  ambien*  noise  and  signal  field. 
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The  method  required  *o  set  the  threshold  level 
•witch**  i*  eumbereom#  and  highly  inaccurate 
since  It  »<■  quire*  amplitude  measurement  of  the 
FUhtr  and  Wiener  output*  on  th*  Develocordar 
di*pUy. 


Recommendation*  for  •  impiilyirig  **tup  and  operation  of  th* 

•  ytt«m  follow.  Th*  thr**hold  level  adjustment  problem  1*  discussed  in  more 
detail  in  Svctioo  IV  and  the  appendix. 


C.  RECOMMENDATIONS 


In  order  to  simplify  initial  eet-up  and  operation  of  the 
Auxiliary  Processor  eystem,  the  following  change*  are  recommended  for 
possible  future  eyeteme.  Moet  of  th*«e  could  be  easily  Incorporated  Into 
the  existing  hardware  with  only  minor  modification#.  Th#  threshold  de¬ 
tector  change  ie  considered  eaaential  for  proper  operation  of  the  eyetem 
as  an  automatic  detection  device.  Recommended  changes  are 


Setting  the  Etcher  transform  constant  Ni  to  a  fixed 


v«iu*  of  111 m  in  order  that  th*  maximum  Fisher 


output  value  wU!  equal  the  maximum  output  number 
(9  bits).  Gain  control  would  aitil  be  available  In  the 
Flaher  D-A  switch. 


Setting  th#  Fisher  transform  constant  Hi  to  a  fixed 
value  of  20,  since,  from  what  is  now  known  about  th* 
expected  sia#  of  F  for  array  data,  this  value  insures 
adequate  number  significance  between  true  Fisher 
values  and  transformed  Either  output  value*  (Figure 

m-i». 


Providing  a  t*-#t  mode  for  displaying  the  Fisher  inter 
mediate  lerr\»  IO,  KI  and  1U  In  order  to  facilitate 
setting  the  Fisber  Input  and  intermedia!*  data  IiskiI 
switches.  The**  terms  could  be  displayed  on  the 
Wiener  power  outputs  whan  in  the  test  mode. 


Prr/tdiag  a*  input  channel  #  sled  vox  capability  for 
the  Flaher  input.  At  present  all  I4CF  input  channel* 
used  U  any  beam-steer  output  are  need  in  the  r  com- 
puts t Son.  For  example,  if  it  war*  desired  to  input  % 

J -component  sensor  to  the  IdCF  for  the  eurpos*  o#  fre¬ 
quency  filiermg,  delaying  and  displaying  tdja<  ent  ii  the 
MCF  deta,  thee*  input*  would  at  pee# be  included  In 
the  F  computation. 
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SECTION  IV 


EVALUATION  OF  ON-LINE  PROCESSING 

The  >rimary  objective  of  the  Auxiliary  Processor  on-line 
evaluation  phase  waa  to  determine  the  feasibility  and  effectiveness  of  on-line 
automatic  detection  processing.  As  an  integral  part  of  this  evaluation,  it 
was  intended  that  data  be  developed  demonstrating  the  interdependence  of  de¬ 
tection  threshold  and  false-alarm  rate.  From  this  data  the  two  detection 
processes,  Fisher  analysis  of  variance  and  Wiener  power,  could  be  compared 
for  detection  capability.  Data  from  the  two  UK  outputs  were  to  be  studied  using 
the  British  Classification  scheme. 

Tho  following  paragraphs  describe  in  detail  the  on-line  evalua¬ 
tion  phase.  For  continuity,  a  description  of  the  processor  operating  configur¬ 
ation  is  presented  first.  Data  for  this  phase  was  developed  primarily  by- 
station  per  onnel  as  a  part  of  the  routine  observatory  operation. 

A.  AUXILIARY  PROCESSOR  PROGRAM 

The  Auxiliary  Processor  was  operational  or  line  at  CPO  from 
30  December  1966  to  17  April  1967.  Table  IV-1  presents  a  summary  of  the 
fixed  program  employed  in  the  digital  MCF  and  Auxiliary  Processor  during 
this  time  period.  Table  IV-2  presents  a  summary  of  the  MCF's  and  beam 
steers  employed  in  the  piocessor. 

As  indicated  in  Table  IV- j,  the  threshold  levels  were  varied 
daily.  Integration  gate  lengths  for  all  outputs  were  initially  programed  for 
2-sec,  but  were  changed  to  3-sec  (60-pts)  on  14  February  1966.  This  change 
in  gate  length  computation  wag  based  upon  a  study  of  mean- signal-length  con¬ 
ducted  in  Dallas.  The  Fisher,  Wiener  and  UK  computations  are  optimised 
when  the  gate  length  of  integration  equals  the  P-wave  duration,  A  mean 
P-wave  duration  of  2.9  sec  was  determined  (SectionV). 


B.  DETECTION  PROCESSING 


Resuks  from  on-line  analysis  of  the  detection  outputs  .ndicate 
that  *he  Fisher  and  Wiener  processes  were  properly  implemented  and  that  this 
type  '■Y  on-line  computation  may  be  highly  effective  in  automatic  de¬ 

tection.  Quantitive  data  demonstrating  the  effectiveness  of  the  present  system 
could  not  be  obtained  due  to  the  effect  of  ron-time  stationary  Fisher  and  Wiener 
Power  noise  distribution  on  the  ttreshold  detection  outpute. 


Table  IV- 1 

PROCESSOR  PROGRAM  DATA 


MCF  PROCESSOR 


Channels 
Filter  pointe 
Multichannel  filter* 
Beam-steers 
Signal  conditioner 
D-A  converter 


19 

57 

5 

9 

All  O’e 

Channel*  0  to  4  *  -3 
Channel*  1 1  to  14  *  -3 


B*am-*teer  hi*tory  15 
Time  delay  28 


AUXILIARY  PROCESSOR 
Arithmetic  drawer 


Output  drawer 


D?»ta  truncation 
•witche* 

Fisher  eummation 
truncation  switch 

Note'.  UK0  and  UK1  are  absolute 
varied  on  a  daily  ba*i** 


UKO 

56 

UK1 

78 

UK  history 

60 

Fisher  N1 

m8 

Fisher  N2 

208 

Fisher  history  60 

MCF  power 

R  *  60,  S  =  1 

UKC 

-3 

UX1 

-3 

MCFO 

-4 

MCF1 

•4 

MCF2 

-3 

MCF3 

-3 

Fisher 

-15 

Fisher 

-3 

UK 

•  2 

MCt.  power 

-2 

-7 


magnitude. 


Threshold  switches  are 


_ _ 


Table  1V-2 

PRESPNT  PROCESSOR  OPER  \T  JNG  h  DDE 


Tide 

Description 

BSO* 

Straight  sum  {21  *  Z  19) 

BS1 

North,  velocity  ■  12.6  km/sec  (Z1  •  Z19) 

BS2 

East,  velocity  *  12.6  km/sec  (Z1  -  Z19) 

BS3 

South,  velocity  ■  12.6  km/ssc  {Zl  -  Z19) 

BS4 

West,  veloctiy  *  12.6  km/ssc  {Zl  -  Z19) 

BS5** 

In-line  summation  toward  USSR  using 

Z4,  6,  7,  9,  14,  15,  17,  and  19 

BS6** 

Transverse  summation  toward  USSR 
perpendicular  to  3S5 

BS7** 

In-line  summation  approximately  toward  NTS 
using  Zl,  2,  8,  10,  11,  12,  13,  and  18 

BS8** 

Transverse  summation  approximately  toward 
NTS  perpendicular  to  BS7 

MCFO 

0.75  cps  low-cut  filter 

MCF1 

MCF3*  convolved  with  1.0-  to  2.  0-cp«  bandpa 
filter 

MCF2 

IFIOWGS*  convolved  with  1.0-  to  2.  0-cps 
bandpass  filter 

MCF3 

MCF3* 

MCF4 

MCF24* 

*  MCFO  subsection  prefilters  the  19-channel  data  and  presents  it  for 
use  in  the  BS  subsection. 

^Used  as  inputs  to  UKO  and  UKl  of  Auxiliary  Processor. 

*  A  complete  description  of  the  filters  is  discussed  in  Appendix  A  of 
CPO  Annual  Report  No.  1. 


1  Automatic  Detection 


Two  ba»ic  problems  whicli  arose  at  CPO  in  attempting  to  imple¬ 
ment  and  evaluate  automatic  on-line  P-wave  detection  with  the  Auxiliary  Pro¬ 
cessor  system  are 

* 

a  The  defined  detection  threshold  for  a  fixed  false-alarm  rate 
was  highly  variable.  Variations  as  great  as  3  to  6  dh  were 
observed  during  a  single  day. 

•  The  method  used  for  rapid  determination  of  the  threshold 
levels  for  automatic  detection  was,  at  best,  inaccurate  and 
cumbersome. 


The  two  problems  actually  interrelate  to  the  extent  that  the  time 
variability  problem  could  have  been  partially  overcome  by  adjusting  the  detec¬ 
tion  threshold  settings  often  (each  hour),  but  a  rapid  and  accurate  technique  for 
determining  the  desired  threshold  level  for  a  fixed  false-alarm  rate  was  not 
available  with  the  existing  hardware  and  setup.  The  inaccuracy  of  the  threshold 
determination  techniques  could  have  been  overcome  by  a  slow  process  of  re¬ 
cording  outputs  on  FM  tape  and  conducting  off-line  analysis  to  determine  the 
desired  levels,  but  this  was  useless  due  to  the  non-time  stationarity  of  the 
Fisher  and  Wiener  power  output  noise  distribution. 

a.  Threshold-Levei  Problem 

A  procedure  was  adopted  at  CPO  to  update  the  threshold  levels  on 
a  daily  basis  (or  more  often,  if  required).  This  procedure,  although  inaccu¬ 
rate,  did  establish  an  approximate  representation  of  the  output  noise  distribu¬ 
tion  in  terms  of  the  daily  threshold-level  setting.  Threshold  levels  for  the 
Fisher  output  and  the  four  Wiener  power  outputs  are  shown  in  Figure  IV- 1  for 
most  of  the  period  the  processor  operated  at  CPO.  The  level  variability  is 
evident  in  this  figure. 

The  exact  cause  of  the  variation  in  output  levels  io  not  known  at 
this  time;  however,  it  is  reasonable  to  suspect  several  sources.  In  the  case 
of  the  Fisher  process,  changes  in  the  mantle  P-wave  noise  level  and  in  the 
level  of  the  trapped-mode  noise  are  probably  the  most  significant  contributors. 
Mantle  P-wave  noise  variability  as  well  as  changes  in  trapped-mode  noise 
direction  could  effect  the  Wiener  detection  outputs. 

Results  from  the  CPO  noise  analysis4  indicated  that  the  predom¬ 
inate  trapped-mode  noise  components  remained  stationary  over  the  analysis 
period.  However,  indications  were  that  the  intensity  of  the  various  compon¬ 
ents  varied;  such  a  variation  *vould  reasonably  affect  the  Fisher  distribution. 


*  As  used  in  this  report,  fixed  false-alarm  rate  is  defined  as  the  number  of 
faise-alarms  per  events  detected  as  a  function  of  time. 
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Recent  work  under  contract  AF  33(657)- 12747  revealed  one 
aource  (the  moat  predominate)  of  m~ntlo  P-wavo  energy  to  be  atmoapheric 
low-pressure  areaa  and  atorm  centera  (hurricanes,  tropical  atorma,  etc.). 

It  ia  reaaonable  to  assume  the  level  of  thia  energy  would  fluctuate  noticeably  as 
a  function  of  distance,  intenaity  and  poaaibly  direction  and  would  aignificantly 
affect  both  detection  output  levels  due  to  the  coincidence  with  aignal  energy  in 
velocity  apace, 

b.  Thre  ahold -Level  Determination  Problem 

Determination  of  the  desired  threahold  level  for  the  Fisher  and 
Wiener  powwr  outputs  would  normally  be  baaed  upon  a  fixed  falae-alarm  rate 
which  could  be  determined  from  knowledge  of  the  detection  output  dia-ribution 
values,  assuming  that  the  distribution  remains  stationary  over  a  reaaonable 
period,  such  as  a  season.  Under  thia  assumption,  off-line  processing  of  de¬ 
tection  output  data  would  be  reasonable  to  determine  the  amplitude  distribution 
and  subsequently  tha  desired  threshold.  Detection  output  data  could  be  collec¬ 
ted  on  FM  tape,  digitized  and  the  distribution  computed. 

Since  the  apparent  amplitude  diitributiona  at  CPO  were  non-time 
stationary  over  perioda  aa  short  as  hours,  ruling  out  off-line  processing  for 
level  determination,  an  on-line  procedure  was  established  which,  in  effect, 
was  a  trade-off  between  accuracy  and  speed.  This  procedure,  although  per¬ 
formed  daily  and  more  often  if  required,  proved  to  be  inadequate  to  handle  the 
non-time  stationanty  problem.  Probebly  thia  due,  in  part,  to  inaccuracy. 

The  procedure  determined  the  mean  amplitude  value  from  a  5- 
mln  period  of  noise  using  five  noise  measurements  at  approximately  1-min  in¬ 
tervals.  Each  measursment  was  made  over  10-see  of  data  and  involved  de¬ 
termining  the  amplitude  value  of  the  peaks  and  troughs  and  arranging  them  to 
obtain  a  mean  output  level.  The  five  measurement!  were  arranged  to  deter¬ 
mine  the  daily  mean  value  for  each  detaction  output.  The  mean  waa  then  mul¬ 
tiplied  by  a  constant  to  determine  the  threahold  value.  Several  constants  were 
tried,  baaed  upon  an  assumed  distribution,  with  the  constant  2,  5  determined 
hast  from  observed  data. 

In  practice,  this  technique  waa  limited  by  the  absence  of  an 
accurate  vertical  scale  on  the  Develocorder.  A  make- shift  scale  waa  applied 
to  the  film  by  placing  the  output  D-A  switches  for  the  Fiaher  and  Wiener  power 
outputs  in  the  M4"  and  "0"  test  modes.  This  produced  a  aero  and  maximum 
value  (777*)  on  the  film  which  could  be  used  ss  a  scale.  When  in  this  test 
mode,  Fisher  and  Wiener  power  data  is  not  output. 

c.  h  ^commendations 

Automatic  detection  could  be  properly  implemented  if  the 
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threshold  lev  ils  coula  be  varied  accurately  and  rapidly  to  compensate  for  die- 
tr.'bution  changea  in  the  Fisher  and  Wiener  power  outputs,  A  suitable  approach 
would  be  the  incorporation  of  an  adaptive  threshold  device  into  the  existing  ard- 
ware  which  could  be  made  to  adapt  to  data  history,  thug  insuring  s  relatively 
consistent  falae-alarm  rate.  By  using  an  on-line  adaptive  algorPhm  to  update 
the  threshold  levels  constantly,  problems  in  measuring  distribution,  etc.  could 
be  avoided.  Only  the  specification  of  a  constant  multiplier  to  be  applied  to  seme 
output  property  such  as  the  mean  value  would  be  required. 

The  appendix  presents  an  adaptivt  thrsshold  algorithm  which  is 
recommended  for  incorporation  into  ths  existing  Auxiliary  Processor, 

2.  Visual  Analysis 

Since  problems  were  encountered  using  the  automatic  detection 
thresholds,  an  attempt  wa»  made  to  analyse  the  Wiener  power  and  Wisher  out¬ 
puts  for  detection  purposes.  This  proved  to  be  unsatisfactory  to  the  station 
analysts  because  the  effect  of  smoothing  of  the  integration  gat*  d*atroy»d  signal 
energy  content  and  waveform,  and  made  the  transition  from  noise  to  signal 
Areas  for  sm*ll  events  quite  gradual  contrasted  to  the  normal  sharp  break  de¬ 
tectable  for  a  p-wave  arrival.  Also,  the  "event  signatures"  which  were  com¬ 
monly  used  to  detect  low-level  signals,  were  not  easily  recognised  after  detec¬ 
tion  processing. 

Xn  addition,  the  absence  of  a  aero  .aval  for  each  of  ths  traces 
made  it  difficult  to  determine  visually  relative  amplitude  on  the  trace  data. 

3.  Comparison  of  Fisher  and  Wiener  Power 

Limited  comparisons  of  the  Fisher  and  Wiener  power  processing 
at  CPO  were  made  by  the  observatory  analysts  while  ettemotirg  to  analyu*  this 
data  for  event  detection,  including 

e  For  high  velocity  signals,  the  ratio  of  the  peak-signal 

output  to  the  RMS  statistic  is  larger  for  the  MCF  processes 
than  for  the  Fisher  output  (Figures  IV-2  and  IV- 3.) 

a  The  Fisher  signel  response  to  quarry  blasts  shows  the 
Fisher  outputs  to  be  unaffected  by  P-  and  energy 

felling  in  the  velocity  ranges  6. 1  to  8.  ktnffMi  and  3.  28 
tf  3.81  km/sec,  respectively  (Figure  IV-4). 

The  Fisher  signal  response  property  was  also  investigated  in 
off-line  research  (Section  V)  and  found  to  decrease  rapidly  with  increasing  wave 
number.  This  property  proved  quite  useful  at  CPO  fo"  distinguishing  local 
and  near  regional  events  from  telessismie  P-wave  energy  end  on  several  oc- 
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Figure  IV-2,  Data  Processed  by  the  MCF  and  Auxiliary  Systems  During 
Known  Signal  Conditions 
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casinos  aid*d  in  d*t**'tlng  t*l***l*mi<  f**wa**  onargf  during  euarry  b«#*t 


irnvaii 


C,  UK  PKOCCS4IHO 


Analyst*  of  (h«  (wq  UK  eu< pvt #  at  CP©  iodic*t*d  ihm  iH» 

Brittah  Classification  #ch* n>*  <e«W  b*  *ff*cti**iy  *bpU*d  oo-l.n*  Ho»*u*r, 
(hr a*  ItmiUiiom  *xiats  two  tay**d  by  th*  pf«*«rt  »»d  iM  thar 

rimed  b*  CP©  non*  prep*rfi*s  Th*»*  UmiUJirsi  *r« 


*  Ths  !«o  UK  out pu! •  ar*  computed  fr  rn  a  ft**d 
progiam  which  In  ***«  t*  r*#tricta  »U**UWMIM 
worii  la  Iwa  aignal  sourc**,  At  CPO  th*  ouJj‘  il *  w*r« 
computed  lor  NTS  and  for  Bo”*!*, 


C'*hiI)c*U‘js  butd  upon  th*  UK  sehstm* 
prsssrvauoe  of  oiff!*!  wiwforifi  Th*  MCF  Auxt* 
t i •  r y  Pro*  a itor  *y*t*rn  I*  limJlod  to  ?l»4i  dynamic 
ran g*  on  lepaf  {11  but).  Sines  Ihl*  ay»i*m  I* 
uttd  priminty  for  dturiloa  bind  upon  ■  ipprasa  on 
of  c oh* real  noli*.  lb*  Input  noli*  l***l  moil  b*  suf- 
fic 1**51)*  high  to  tnaur*  adequat*  noli*  alatiatle*  for 
MCF  proc***ing.  Thu*.  I*rg*  signals  *r*  clipped 
on  input  or  during  iW* rrr *diat*  computations. 


Tb*  UK  t*th»»^«*  I*  designed  for  application  to  |«. 
t«rm*4i*t*  **a*  ereaaarray  data  (appro*  ima  isrly 
10  km)  to  in*ur*  ad*quat*  i* solution  and  an  oncer 
r#lat*4  ambi*it  nei**  fi«W.  Th*  CFO  array  {S,t 
km  dli)  violated  th*  ittiunpIlM  of  oncer  elat»« 
nou*  and  provid*d  vary  poor  directional*  velocity 
resolution 


At  CP©  th*  UK  portion  of  th*  prwitwr  **i  pro|r»fn*d  for 
mammon  *t«<orm  pr***f*ill«i  within  th*  limitation  of  b*A*|  th*  input 
dis*  gain  salting  optimised  for  MCF  programing  Thus  th*  UK  input  data 
truncat. on  ***itbei  »*r*  set  such  that  a  li*bu  aignal  an  input  *e  the  MCF 
{nuuumom  P*P  wUhevt  clippmgj  would  b*  processed  by  th*  UK  .artipjUUon 
and  output  without  «l»pp  ng.  Thl#  s*i.«p  *u  considered  optimum  for  th* 
tUiiifittiiM  cempuiaJioii  tine*  oivtftrj.*  preservation  I*  men  ImpcrUM. 
Howwv*r«  with  th*  D***locer4*r  adjusted  to  a  measurable  gain  .aval,  moduln- 
tion  '  f  fmv4***i  aignal  *»*  poor. 


Examples  of  UK  output  data  for  th*  two  programed  signal  re- 
gi-  r  •  in  shown  in  Figures  IV  •  V  through  IV*  -  these  figure*  demonstrate 
ih*  dynan  »c  rang*  limitation  Figure  |¥*7  presents  an  example  s!  a  quarry 
blast  processed  by  the  UK  technique,  For  this  blast  and  for  the  special  event 
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SECTION  V 


SUPPORTING  DETECTION  PROCESSING  RESEARCH 

In  conjunction  with  on-line  operation  of  the  Auxiliary  R-ocessot, 
off-line  reaearch  wit  dt  ne  in  Dallas.  Primary  reaearch  goala  were  to  verify 
the  logic  of  the  Auxiliary  Proceaaor,  to  determine  and  atudy  and  parametera  usod 
with  the  different  proceaaea,  and  to  atudy  the  effecta  cf  varying  theae  para¬ 
meter. 

A.  PROCESSOR  SIMULATION  PROGRAM 

A  aimulation  program  to  sf.uJ’  the  proceaaor  logic  waa 
written  which  used  the  IBM  7044  computer.  Thia  program  aimulated  the  logic 
of  >  th  the  MCF  and  Auxiliary  Processor  and  checked  for  round-off  errors  in 
the  proceaaor  computations.  A  complete  description  of  the  program  may  be 
found  in  CPO  Special  Report  3 

B.  THEORETICAL  NOISE  SAMPLE 


To  provide  i  epreaentative  CPO  ambien..  noise  statistics  for 
use  in  studying  the  Fisher  statistic  while  minimizing  computer  requirements 
(i.  e.,  number  of  noiae  samples  required  for  processing)  a  5000-pt  theoretical 
noise  sample  waa  synthesized.  This  sample  was  designed  to  have  the  same 
correlation  matrix  as  the  CPO  noise-ensemble  samples  A,  B,  E,  F,  and  I, 
developed  under  Contract  AF  33(657)- 12331.  Due  to  the  noiae  field  time 
stationarity  (CPO  Special  Report  No.  I4),  it  was  felt  that  this  1963  average 
noiae  ensemble  would  provide  a  statistical  representation  of  the  CPO  ambient 
noise  field. 

The  technique  used  to  develop  this  noise  sample  (Figur^  V-l) 
is  discussed  in  Section  III  of  the  Array  Research  Semiannual  Report  3  .  After 
the  data  were  developed,  the  noise  sample  was  tested  for  the  same  properties 
as  the  desired  correlation  matrix.  Power  spectra  were  computed  from  random 
correlations  using  both  tht  theoretical  and  original  noise  sani'ples.  These 
power  spectra  (Figure  V-2>  show  that  the  theoretical  noise  sample  has  the 
si  me  properties  as  the  original  noise  sample  and  does  represent  true  seismic 
data. 

C.  PROCESSING  PARAMETERS 

Choice  of  correct  operating  parameters  is  of  primary  importance  in  on¬ 
line  operation  of  the  Auxiliary  Processor.  Research  in  Dallas  has  resulted  in 
optimizing  two  critical  parameters:  the  signal  gate  length  and  corner  frequency 
of  MCFO  which  is  the  prefilter  for  the  Fisher  process. 
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Figure  V-i.  CPO  Theoretical  Noiie  Sample 
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I.  Gate  Length 


To  establish  tha  optimum  gate  length  for  the  Fisher  and  Wiener 
compulations,  a  study  was  made  to  determine  the  time  duration  of  the  primary 
P-wave  signal  from  teleeeisms  recorded  over  a  6-day  period.  The  processes 
are  optimised  when  the  computation  gate  length  equals  the  signal  duration. 

Data  for  this  study  came  from  a  visual  interpretation  of  the  CPO  Devnlocorder 
film.  P-wave  duration  from  57  teleseismic  signals  was  measured,  yielding  a 
mean  P-pulse  signal  duration  of  2.90  sec  and  a  variance  In  P-pulse  signal  du¬ 
ration  of  1 . 90  sec. 

Upon  completion  of  this  study,  the  Auxiliary  Processor  program 
for  the  Fisher,  Wiener  power  and  UK  processas  was  changed  from  2.0  to  3.0 
sec.  After  this  adjustment,  analysts  reported  data  easier  to  analyse,  since  the 
processed  traces  (notkbly  ths  Fisher)  were  stabilised  and  did  not  correlate  as 
well  in  time. 

2.  Optimum  Low-Cut  Filter 

Several  low-cut  filtere  with  varying  corner  frequencies  were 
developed  to  determine  the  optimum  low-cut  filter  to  be  used  in  the  MCFO  sub¬ 
section  of  the  processor.  The  choice  of  this  filter  is  extremely  important 
since  it  is  the  prefilter  for  the  Fieher  subroutine.  Filtering  of  the  input  data 
to  the  Fieher  process  is  necessary  in  order  to  remove  highly-correlated  low- 
frequency  microseismic  energy  so  that  the  Fisher  assumption  of  spatially  un¬ 
correlated  energy  is  approximated. 

Analysis  of  filters  applied  on-line  at  CPO  (Figure  V-3)  showed 
the  optimum  corner  frequency  to  be  approximately  1.0  cps.  The  filters  with 
0. 75-,  1.0-,  and  1,25-cpt  corner  frequencies  were  subsequently  applied  off¬ 
line  to  the  theoretical  noise  sample  (Figure  V-l)  and  the  Fisher  output  com¬ 
puted  for  each  of  the  three  case*.  Hesults  of  this  processing  are  shown  in 
Figura  V-4,  To  facilitate  the  interpretation  of  this  data,  a  cumulative-distri¬ 
bution  counter  wac  added  to  the  Fisher  subsection  of  the  simulation  program. 
Th"*  cumulative  distributions  for  this  study  are  shown  In  Figure  V-5. 


•When  interpreting  Figure  V-5  and  the  other  cumulative  distributions  presented 
in  this  report,  the  horisontal  scales  are  not  linear.  Even  though  the  eurree  may 
appear  very  similar,  their  true  meaning  can  be  interpreted  by  noting  the  last 
saveral  inchae  on  tha  horlaentai  scale. 
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Figure*  V-4  and  V-5  »how  that  the  optimum  low-cut  prefilter  *hould  have  a 
corner  frequency  of  0.  75-1. 00. 


G 'eater  suppression  of  coherent  noi«e  could  be  obtained  by 
choosing  a  filter  with  a  corner  frequency  higher  than  0.  75  cps,  but  this  would 
emit  extensive  signal  degradation  since  the  predominate  portion  of  the  F-wava 
usually  occurs  between  0.6  and  1.0  cps. 

The  average  signal  spectrum  was  investigated  to  determine  the 
optimum  trade  off  between  minimizing  the  Fisher  noioe  distribution  (low-cut 
frequency  filtering)  and  maintaining  adequate  s*  energy.  The  average  CPO 
signal  spectrum  was  computed  for  an  ensemble  ol  Kurile  Islands  events  de¬ 
veloped  under  Contract  AF  33(657)- 12747  8  and  is  shown  in  Figure  V-6. 

Peak  signal  energy  occurs  at  0.  68  cps.  since  this  result  could  be  tuned  to  the 
Kurile  Island  ■*  region,  and  not  represent  an  average  of  the  CPO  predominant 
P-wave  energy,  additional  P-wav  •  data  was  derived  from  the  CPO  standard 
analysis  forms.  An  average  predominate  P-wave  frequency  of  0.82  cps  was 
determined  from  one  month  of  standard  station  report  data. 

Based  on  the  Fisher  noise  output  distribution  shown  in  Figure 
V-5,  and  with  the  knowledge  that  the  peak-signal  spectrum  lies  between  0.  68 
and  O.  82  cps,  the  0.75  cps  low-cut  filter  (Figure  V-3)  was  chosen  for  on-line 
application. 

D.  PROCESSED  NOISE  SAMPLES 

1.  Threshold  Non- Time  Stationarity 

The  Fisher  output  noite  distribution  non-time  stationarity  was 
studied.  Four  noise  samples,  I,  II,  III,  and  IV,  were  chosen  for  processing 
from  the  1965  CPO  library  (Figure  V-7)  and  represent  the  different  types  of 
noise  backgrounds  encountered  at  CPO  (low-I  medium -II  and  III  and 
high— IV).  Power  deisity  spectra  for  these  samples  are  shown  in  Figure  V  -8. 
The  Fisher  statistic  for  the  four  samples  is  shown  in  Figure  V-9,  and  the 
cumulative  distributions  for  the  F4sh»r  statistics  are  shown  in  Figure  V-10. 
These  illustrations  show  the  Fisher  statistic  to  be  larger  for  high  noise 
(Sample  IV)  and  smaller  for  the  other  samples. 

Although  Sample  I  is  of  low-level  noise,  its  Fisher  statistic  is 
*arger  than  those  from  Samples  II  and  in.  While  the  cause  for  this  is  unknown, 
it  is  possible  that  the  mantle  P-wave  noise  level  was  much  highc-  for  Sample  I 

*Thc  theoretical  noise  sample  (Figure  V-l)  has  been  prewhitened  with  a  1.0 
cps  12  db/octave  low-cut  filter  which  closely  approximates  the  frequency  fil¬ 
tering  accomplished  on  the  MCF  analog  output  data  in  the  MCF  signal  condi¬ 
tioning  lrawer. 
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Figure  V-7.  Hoi.e  Sample.  1,  11,  ill,  and  IV  from  1965  CPO  Noi.e  Library 
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than  for  SamOea  II  and  IK,  caualng  thia  raault. 

Figurea  V-9  and  V-10  ahow  that  the  large  majority  of  Fiaher 
valuea  (the  50-percent  level)  occur  at  approximately  F«2.0.  This  value  cloaely 
agreea  with  that  publiahed  in  earlier  work  on  the  Fiaher  atatiatic.  The  real 
aignificant  of  thia  data  ia  that  the  cumulative  diatribution  in  Figure  V-10  ahowa 
that  the  threahold  level  for  fixed  fulae-alarm  rate  ia  variable  and  not  atrongly 
related  to  the  power-denaity  apectrum  of  the  ambient  noiae. 

2.  Effect  of  Correlated  Noiee  on  the  Fiaher  Statietic 

To  atudy  the  effect  cl  correlated  noiae  on  the  Fiaher  atatiatic 
(which  aeaumea  a  apatially  random  noiae  diatribution),  a  aimple  recorded  in 
1963^  waa  proceaaed  along  with  the  data  preaentod  previoualy  in  Figurea  V-4 
and  V-5.  The  noiae  aample  ia  ahown  in  Figure  V-ll,  and  ita  Fiaher  atatiatic 
ia  plotted  in  Figure  V-12.  Tha  cumulative  diatribution  for  the  aample  la  ahown 
in  Figure  V-13,  uaing  no  filter  and  for  corner  frequenclea  of  0.75,  1.0  and  1.25 
cpa.  Reaulta  are  aimilar  to  thoae  preaented  in  Figure  V-5  and  indicate  th  it 
the  correlated  noiae  atrongly  effecta  the  Fiaher  output  atatiatic  below  1.0  cpa. 
The  amall  difference  in  the  no  filter  and  0.  75  cpa  filter  diatributiona  in 
Figurea  V-5  and  V-13  are  due  to  the  fact  that  the  aample  in  Figure  V-ll  waa 
not  prewhitened. 

Attempta  to  correlate  the  diatribution  changea  to  noiae  pre¬ 
dictability  were  unaucceaaful,  Preaent  eatimatea  of  the  mantle  P-wave  noiae 
■  peetrum  indicate  that  thia  energy  ia  predominate  at  lower  frequenciea  and 
rapidly  decaya  up  to  1.0  -  1.4  cpa  which  could  aignifieantly  affect  the  Fiaher 
and  Wiener  outputa  below  1.0  cpa. 

E.  SIGNAL  PROCESSING 

A  Bet  of  theoretical  white  aignal  waveieta  waa  generated  to 
determine  the  relation  between  Fiaher  output  value  and  aignal  velocity.  Theae 
waveieta  were  generated  for  infinte  apparent  horizontal  velocity  ana  apparent 
horizontal  velocitiea  of  25,  12  and  8  km/aec.  Azimuthal  directiono  were 
choaen  ao  that  the  aignala  would  be  incident  to  the  array  at  pointa  of  maximum 
(180°)  and  minimum  (27U°)  valuea  of  the  array  reaponae.  The  waveieta  were 
aynthetically  added  to  the  theoretical  noiae  aample  (Figure  V-l)  and  the  Fiaher 
ititiitic  computed.  Figure  V-14  ia  an  example  of  the  noiae  p*ua  wavelet  for 
the  infinite  velocity  caae.  Due  to  the  aimilaritiea  between  the  plota,  only  the 
cumulative  Fiaher  diatributiona  are  ahown  in  Figure  V-15,  From  Figure  V-15, 
two  concluaiona  are  drawn. 

e  For  a  given  velocity,  the  direction  haa  very  little  effect 
on  the  Fiaher  valuea.  Thia  reault  ia  expected  due  to  the 
near  aymmetry  of  the  CPO  array. 
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•  The  size  of  the  Fisher  values  rapidly  increase  as  Me  velocity 
of  the  incident  signal  increases. 

Tne  decrease  in  signal  output  level  for  the  Fisher  computation 
is  greater  chan  would  be  expected  on  the  basis  an  array  response  comparison, 
and  explains  the  Fisher  property  of  significantly  suppressing  quarry  blast  in¬ 
formation  (Section  III).  When  compared  with  Wiener  pwoer  processing,  the 
signal  suppression  property  is  more  severe  for  the  Fisher  output.  The  Wiener 
power  output  response  is  easily  determined  from  the  filter  response,  which 
for  the  CPO  MC/'s  in  Table  III-  1  at  1.0  cps  is 


Infinite 


Filter 

Velocity 

25  km/ sec 

12  km/sec 

8  km/sec 

MCF1 

-0  db 

-1  db 

-4  db 

-5  db 

MCF2 

+3  db 

+2  db 

-0  db 

-4  db 

MCF3 

-0  db 

-1  db 

-4  db 

-5  db 

MCF4 

-3  db 

-3  db 

-3  db 

-4  db 

Maximum  response  variation  in  the  MCF  data  is  7  db  for  MCF2, 
while  the  Fisher  varies  12  db  over  the  same  velocity  range  (infinite  velocity 
to  8  km/sec). 


Fisher  output  signal  degradation  is  partially  implied  from  th< 
observation  made  in  Section  HI  regarding  the  peak-signal-to-RMS-noise  output 
comparison  with  the  MCF  processes, 

Fisher  output  responre  can  be  maximized  in  practice  if  the  input 
data  is  beam-steered  to  enhance  signal  energy  prior  to  computation  of  the 
Fisher  statistic.  This  approach,  although  well  suited  to  off-line  processing, 
would  be  impractical  with  the  existing  hardware  because  of  the  fixed  program 
restriction. 
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APPENDIX 


ADAPTIVE  THRESHOLD  DETECTION 

At  present,  use  of  the  existing  MCF  Auxiliary  Processor 
system  for  automatic  detection  is  limited  due  to  non -time  stationary  Fisher 
and  Wiener  power  output  noise  distributions.  This  1  nitation  could  be  over¬ 
come  by  incorporating  adaptive  threshold  detection  into  the  Auxiliary  Pro¬ 
cessor. 


An  adaptive  threshold  detector  is  one  which  automatically  ad¬ 
justs  the  detection  threshold  levels  based  upon  the  output  noise  history  as  de¬ 
fined  over  some  specified  period.  Associated  with  the  automatic  detection 
capability  must  be  a  specified  false-alarm  rate  and  a  learning  rate. 

Since  a  modification  to  the  existing  hardware  is  recommended, 
the  storage  requirements  for  the  adaptive  algorithm  should  be  minimized,  and 
the  modifications  should  be  incorporated  into  the  present  Auxiliary  Processor 
drawer  space. 


A.  RECOMMENDED  SOLUTION 


In  order  to  meet  the  assumed  hardware  restrictions  while  still 
providing  an  adequate  technical  solution,  it  is  recommended  that  the  threshold 
level  be  determined  from  output  history  data  which  has  been  smoothed  by  an 
exponential  time  window.  Mathematically,  the  output  Fisher  or  Wiener  power 
level  would  be  described  by  a  function  X j  +  j  given  ^y: 
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whe*  e: 


i  *  time 

j  =  output  channel  (one  for  Fisher  and  four 
for  the  Wiener  power) 

X'1  s  Outnut  value  of  the  jth  channel  at  time 

i-n+1  , 


time  constant 
characteristic 


This  system  provides  for  control  of  the  learning  speed  with  the 
a,  and  includes  all  past  history  data  to  determine  the  ©utput 
.  The  threshold  would  then  be  given  by  some  value 
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where  K  is  a  constant  which  sets  the  desired  falae-alarm  rate. 


The  exponential  smoothing  technique  provides  for  weighting  of 
the  output  statistic  x|  +  j  as  a  function  of  time,  i.e.,the  most  recent  data  is 
weighted  heaviest. 

B.  IMPLEMENTATION 


follows: 


The  above  equation  can  be  simplified  for  on-time  application  as 


letting 


_  j  r 

X  .  ,  =  L  A  XJ  e 

1  +  1  n=0  i - n  +  1 


-na 


=  X' 


J  .*L 


-na 


i  +  1  n  =  1  i -n  +  1 


n  =  n  +1 


x.jA .  =  X^  +  Z.  X,  , 

*  +  1  i-fl  n-0  i-n 


-(n'  +  l)a 


*  X 


i  +  1 


-a  7  j 
te  £'*0Xi-n' 


■na 


Howeve 


r,  xj  *  I  x|  ,  e"n’a 

i  n'  *  0  i  -  n' 


v  j  •  vi 


Therefore,  XJ  .  »  XJ  .  +  X 
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and  the  updated  threshold  is  given  by 
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The  threshold  computation  for  the  one  output  is  thas  reduced  to 
two  multiples  and  one  add  operation  and  storage  requirements  are  reduced  to  a 
single  word,  X  j  . 

A  preliminary  investigation  of  the  hardware  requirements  for 
implementing  this  adaptive  threshold  algorithm  indicates  that  existing  hard¬ 
ware  can  be  modified  to  include  this  capability  with  no  increased  core  storage 
or  drawer  space  requirements. 
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